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Glacial cycles of the Quaternary have heavily influenced the demographic history of various species. To test 
the evolutionary impact of palaeo-geologic and climatic events on the demographic history of marine taxa 
from the coastal Western Pacific, we investigated the population structure and demographic history of two 
economically important fish ( Trichiurus japonicus and T. nanhaiensis) that inhabit the continental shelves 
of the East China and northern South China Seas using the mitochondrial cytochrome b sequences and 
Bayesian Skyline Plot analyses. A molecular rate of 2.03% per million years, calibrated to the earliest 
flooding of the East China Sea shelf (70-140 kya), revealed a strong correlation between population sizes 
and primary production. Furthermore, comparison of the demographic history of T. japonicus populations 
from the East China and South China Seas provided evidence of the postglacial development of the 
Changjiang (Yangtze River) Delta. In the South China Sea, interspecific comparisons between T. japonicus 
and T. nanhaiensis indicated possible evolutionary responses to changes in palaeo-productivity that were 
influenced by East Asian winter monsoons. This study not only provides insight into the demographic 
history of cutlassfish but also reveals potential clues regarding the historic productivity and regional 
oceanographic conditions of the Western Pacific marginal seas. 



Glacial cycles and the consequent sea-level fluctuations of the Quaternary epoch greatly influenced the 
spatial distribution and population dynamics of a number of terrestrial 13 , freshwater 4,5 and marine 
species 6 ". This was particularly true for marine organisms of oceans and marginal seas 911 . The advan- 
cing terrestrial ice sheets and falling sea level during glacial periods exposed shallow coastal continental 
shelves 1214 , greatly reducing habitat availability for many marine taxa. However, rising sea levels during the 
interglacial periods inundated continental shelves and enlarged marginal seas on a global scale 15 . These cyclic 
changes in habitat availability likely led to repeated range expansions and contractions, influencing the popu- 
lation dynamics of various taxa 1617 . 

The China seas, including the South China Sea (SCS), East China Sea (ECS), Yellow Sea (YS) and Bohai Sea 
(BH), are relatively young (i.e., <30 million years [myr]) marginal seas of the Western Pacific Ocean 18 (Fig. 1). 
The Bohai, Yellow and East China Seas share a wide and shallow continental shelf, whereas the South China Sea 
forms a deep sea basin with a narrow continental shelf in the north and the wide Sunda Shelf in the southwest. 
These marginal seas are supplied with high levels of terrigenous nutrients from the Asian continent through river 
run-off, which has led to high primary productivity and increased fish production 19 . During the Pleistocene, sea- 
level fluctuations would have repeatedly exposed and inundated these shallow marginal seas, dramatically 
changing nutrient input and habitat availability. For example, some recent studies have revealed similar patterns 
of population expansions for a number of marine taxa in response to Pleistocene sea-level fluctuations in the 
China seas 10,20 22 . Furthermore, the Yellow-Bohai, East China and South China Seas have been identified as large 
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Figure 1 | Map of the East China and northern South China Seas from 
which individuals of Trichiutus japonicus and T. nanhaiensis were 
sampled for this study. Individual sampling sites are listed in Table 1. The 
shoreline of the LGM ( — 130 m) is indicated by thin dashed line and the 
boundary of the East China Sea and South China Sea is indicated by thick 
dashed line. BH = Bohai Sea, YS = Yellow Sea, ECS = East China Sea, SCS 
= South China Sea, OT = Okinawa Trough, DS = Dongshan, GT = Gulf 
of Tonkin. Map was created in software Generic Mapping Tools (GMT 
v3. 1. 1; http://gmt.soest.hawaii.edu/) . 

marine ecosystems (LMEs) that have global importance, based on 
their distinct bathymetry, hydrography, productivity and trophic 
relationships 23 , and the different oceanographic conditions of these 
LMEs are expected to have influenced the evolutionary histories of 
regional populations. However, comparative studies investigating 
the dynamic responses of different marine taxa to endemic oceano- 
graphic conditions among these LMEs are lacking. Influenced by 
past geologic and climatic events, understanding the demographic 
history of regional populations of marine organisms from the China 
seas may provide a framework for elucidating the historic climatic 
and oceanographic conditions of the Western Pacific Ocean. 

The Trichiuridae (cutlassfish or hairtails) are an economically 
important group of benthic fish that inhabit the continental shelves 
and shallow seas (20-200 m) of tropical and temperate waters 24 , with 
over 80% of global cutlassfish catch production originating from the 
China seas 25 . Furthermore, cutlassfish are the most important com- 
mercial marine fish species in East Asia, accounting for 8.5-18.3% of 
the total annual fishing catch in the China seas between 1950 and 
2000 26 . Given these pressures, it is not surprising that the over- 
exploitation of cutlassfish in recent years has led to a dramatic 
decrease in catch production and population size 27 . Additionally, 
increased fishing pressures have led to reductions in average body 
size and increases in maturation rates 28 . Given the glacial cycles 
mentioned above and these recent anthropogenic pressures, cutlass- 
fish in the Western Pacific region are predicted to have a complex 
and dynamic demographic history. 

The demographic history of populations can often be revealed by 
analysing the genetic variation using coalescent-based statistical 
methodologies, such as mismatch distribution 29 and Bayesian 



Skyline Plot (BSP) analyses 30 . For example, these types of analyses 
have revealed repeated periods of population expansions and con- 
tractions due to glacial cycles, sea-level fluctuations and anthro- 
pogenic pressures for a number of marine taxa such as mitten 
crabs {Eriocheir sensu stricto) in the East Asian marginal seas 31 , 
North Pacific herring (Clupea pallasitf and Antarctic fur seals 
(Arctocephalus gazelle) 32 . In this study, the genetic variation and 
demographic history of two cutlassfish species, T. japonicus from 
the northern South China Sea and East China Sea and T. nanhaiensis 
from the northern South China Sea, were investigated via sequence 
analyses of the mitochondrial cytochrome b (Cyt b) gene and via BSP. 
Specifically, we attempted to investigate the evolutionary influences 
of past climatic events and regional oceanographic conditions on the 
demographic histories of marine organisms, including T. japonicus 
and T. nanhaiensis, in two West Pacific marginal seas, the northern 
South China Sea and the East China Sea. 

Results 

Population genetic analyses. A total of 761 cutlassfish were collected 
from 26 sampling sites in the East China Sea and northern South 
China Sea (Fig. 1). Direct sequencing of the Cyt b gene resulted in 509 
and 252 sequences belonging to T. japonicus (400, East China Sea; 
109, South China Sea) and T. nanhaiensis (South China Sea), respec- 
tively. From these, 239 T. japonicus and 74 T. nanhaiensis unique Cyt 
b haplotypes were identified and deposited into GenBank under the 
accession numbers KF551253-KF551565. Specific sample sizes and 
numbers of haplotypes recovered for each species from each locality 
are summarised in Tables 1 and 2. Estimates of haplotype diversity 
(h) and nucleotide diversity (n) for T. japonicus were mostly consis- 
tent among sampling sites and geographic populations (Table 1). 
However, high haplotype diversities and low nucleotide diversities 
were observed for T. nanhaiensis at each locality, indicating recent 
population expansions (Table 1). 

A star-like network radiating from a central haplotype was recov- 
ered from the median-joining network analyses of T. japonicus Cyt b 
sequences (Fig. 2a). Within the network, haplotypes sampled from 
the East China Sea and northern South China Sea were evenly dis- 
tributed. Network analyses of T. nanhaiensis Cyt b sequences from 
the northern South China Sea produced one network with two diver- 
gent clusters (I and II) distinguished by nine substitutions (net dis- 
tance, 0.989%; Fig. 2b). The majority (243) of T. nanhaiensis Cyt b 
sequences belonged to cluster I and exhibited a diffuse network pat- 
tern deriving from three dominant haplotypes. The remaining nine 
sequences belonged to cluster II. No obvious geographic structure 
was observed between clusters I and II. 

Molecular rate calibration. Using the equation |i = i/(2mt) and the 
flooding of the East China Sea shelf 70- 140 kya as a calibration point, 
a mutation rate of 2.03% (1.35-2.7%) per myr was inferred. If the 
earliest population expansion of cutlassfish in the East China Sea 
induced by the postglacial high sea levels of <22 kya was used as a 
calibration point, a substitution rate of over 9.7% per myr could be 
inferred, which is approaching the maximum substitution rate 
estimated for the mitochondrial control region of fish 33 . Because 
the substitution rate of Cyt b is much slower than the mitochon- 
drial control region in fishes 34 , the calibration point of the Last 
Interglacial sea level highstand (70-140 kya) used here appears to 
be more appropriate. However, given the uncertainty of rate 
estimates, the maximum {i.e., 2.7% per myr) and minimum (i.e., 
1.35% per myr) estimated rates (bracketed values) were also used 
in the BSP inference analyses below. 

Population dynamics. The pairwise mismatch distributions for T. 
japonicus as a whole and the two geographic groups of T. japonicus 
were unimodal (Fig. 3a, b). Additionally, no significant deviation was 
found between any of the observed distributions and expected 
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1 Sampling locations and indices of genetic diversity of T. japonicus and T. nanhaiensis collected from the South 
Letters indicate variables: n, sample size; nh, haplotype number; h, haplotype diversity; n, nucleotide diversity 
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distributions of haplotypes based on the sudden expansion model 
(Table 2), indicating recent population expansions. Furthermore, 
neutrality tests (Tajima's D and Fu's F s ) revealed significant 
negative values for all species and geographic groups (Table 2). 
Taken together, these results indicate recent population expansion 
in both T. japonicus and T. nanhaiensis and in the two geographic 
groups of T. japonicus. 

Based on extended BSP analyses, the number of groups of coales- 
cent intervals was inferred from each data set with a mean of 2.436 
with 95% credible set {1, 4} for T. japonicus, 2.051{1, 4} for the East 
China Sea population of T. japonicus, 1.598{1, 3} for the South China 
Sea population of T. japonicus and 1.612 {1, 3} for T. nanhaiensis. 
These statistics strongly support a complex population demographic 
model (e.g., BSP) over a single constant population size model (1- 
group model). 

Population size fluctuations through time were further compared 
using BSP among species and geographic groups (Figs. 4, 5), and the 
partitioning of Cyt b sequences between 3 codon positions neither 
improved the resolution nor changed the results (data not shown). 



For T. japonicus as a whole, BSP revealed two episodes of population 
growths (Fig. 4). The first episode of rapid population growth 
occurred 45-105 kya [30-70 kya, 2.7% per myr; 60-140 kya, 
1.35% per myr] and was followed by a reduction in population size 
22.5-37.5 kya [15-25 kya, 2.7% per myr; 30-50 kya, 1.35% per myr] 
that coincided with the Last Glaciation Maximum (LGM). Following 
the LGM and continuing to the present, T. japonicus experienced a 
second and more rapid episode of population growth 0-22.5 kya [0- 
15 kya, 2.7% per myr; 0-30 kya, 1.35% per myr]. The current popu- 
lation size of T. japonicus appears to be two orders of magnitude 
larger than it was at the time of coalescence 304 (202.7-405) kya 
(Fig. 4). 

Despite showing little genetic differentiation or population struc- 
ture, T. japonicus sampled from the East China Sea and northern 
South China Sea exhibit distinct demographic histories. Similar to 
the species as a whole, BSP analysis revealed two episodes of popu- 
lation growth for T. japonicus from the East China Sea (Fig. 4). 
However, T. japonicus from the northern South China Sea only 
exhibited a single period of population growth occurring 45- 



Table 2 Results of dynamic tests and molecular clock analyses for T. japonicus and T. nanhaiensis from the South and East China Seas. 
Letters indicate variables and abbreviations: n, sample size; nh, haplotype number; h, haplotype diversity; it, nucleotide diversity; SSD = 
sum of square deviations between the observed and expected mismatch, t = mutational time scale. *P < 0.05, ** P < 0.001 


Species Group n nh h 


n (%) Tajima's D 


Fu's F s 


Coalescent time 
SSD/P t (kyr, 2.03% per myr) 


T.japonicus Total 509 239 0.9741 
South China Sea 109 73 0.9849 
East China Sea 400 193 0.9726 

7. nanhaiensis Total 252 74 0.8727 


0.55 -2.3359** 
0.57 -2.1519** 
0.54 -2.331** 
0.23 -2.3928** 


-24.4302** 
-25.0221** 
-24.5337** 
-26.3877** 


0.0007/0.614 4.344 303.98(178.86-455.07) 
0.0017/0.541 4.984 265.53(147.20-398.72) 
0.0007/0.659 4.289 299.22(174-452.60) 
0.0021/0.07 1.930 291.11(145.76-457.52) 
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Figure 2 | Evolutionary networks depicting relationships among Cyt b haplotypes recovered from T. japonicus and T. nanhaiensis. Geographic 
information for haplotypes is distinguished using white (East China Sea) and black (South China Sea). The sizes of circles are proportional to the 
frequency at which specific haplotypes were recovered. Numbers indicate mutational differences of two or more substitutions. For each network, red dots 
represent unsampled (missing) haplotypes. 



105 kya [30-70 kya, 2.7% per myr; 60-140 kya, 1.35% per myr], 
similar to that observed in the East China Sea (Fig. 4). This period 
of population growth (45-105 kya) resulted in a 10-fold increase in 
population size for both the East China Sea and northern South 
China Sea populations relative to those at the times of coalescence 
for these populations. However, in the northern South China Sea, 
this growth continued, peaked during the LGM (22-30 kya) and has 
since undergone a steady decline (Fig. 4, Fig. 5). As a result, the 
present population size of T. japonicus in the East China Sea is ca. 
20 times larger than that in the northern South China Sea (Fig. 4). 
Haplotypes of T. japonicus from the northern South China Sea 
exhibit a younger coalescent time of 266 kya [177-354 kya] 



compared with that of 299.3 kya [199.5-399 kya] for haplotypes of 
T. japonicus from the East China Sea. 

T. nanhaiensis from the northern South China Sea appears to have 
responded differently to the climatic dynamics of the area (Fig. 5). 
BSP revealed a steady population decline 37.5-165 kya [25-110 kya, 
2.7% per myr; 50-220 kya, 1.35% per myr] followed by a rapid 
population expansion (15-37.5) [10-25 kya, 2.7% per myr; 20- 
50 kya, 1.35% per myr), with an estimated 100-fold increase in the 
maximum effective population size. However, beginning —15 kya 
[10-20 kya] and continuing to the present, T. nanhaiensis has been 
undergoing a second bottleneck event, and its population has been 
reduced to a size similar to that of the northern South China Sea 
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Figure 3 | Mismatch distribution analyses of T. japonicus and T. nanhaiensis in the China seas, (a) T. japonicus and T. japonicus population of the East 
China Sea; (b) T. japonicus population of the South China Sea and (c) T. nanhaiensis population of the South China Sea. 
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population of T. japonicus (Fig. 5). Haplotypes of T. nanhaiensis in 
the northern South China Sea were estimated to have coalesced 
291 kya [194.1-388.2 kya]. 

Discussion 

Because substitution rates of the mitochondrial Cyt b gene in fish 
appear to be variable 35 39 and time-dependent 40 ""**, the present study 
used two recent historical events to calibrate a molecular clock for 
Trichiurus (Supplementary information SI). Despite some uncer- 
tainty, the molecular rate of 2.03 (1.35-2.7)% per myr estimated here 
provides an approximate timescale for better understanding the 
influence of late Pleistocene climatic events on the evolution and 
ecology of marine organisms. 

As top predators, the population abundance of marine fish has 
been considered a good indicator of primary productivity 45,46 , which 
appears to hold true for the population demographic history of T. 
japonicus and T. nanhaiensis in the China seas. The differences 
between present-day (t = 0) estimated effective population sizes 
for intraspecific geographic groups (East China Sea vs. northern 
South China Sea) of T. japonicus or interspecific geographic groups 
(the East China Sea population of T. japonicus vs. the northern South 
China Sea population of T. nanhaiensis) reach an order of magnitude 
(Figs. 4, 5), suggesting that the East China Sea supports ca. 20 times 
more cutlassfish than the northern South China Sea. These differ- 
ences in the estimated regional effective population sizes are consist- 
ent with cutlassfish surveys in the China seas, where over 80% of the 
annual catch of cutlassfish originates from the East China Sea 26 . 

Several factors may have led to the higher effective population size 
and biomass of T. japonicus in the East China Sea than in the north- 
ern South China Sea. First, as mentioned above, cutlassfish are ben- 
thic fishes that inhabit the photic zone above continental shelves at a 



depth of 20-200 m 24,26 . In this context, the East China Sea has a much 
wider continental shelf (540,000 km 2 ), providing a greater area of 
suitable habitat than the narrower and deeper northern slope of the 
South China Sea (370,000 km 2 ) 47 . Second, the East China Sea pos- 
sesses a higher average annual primary productivity (609.65 gC 
m~ J a _1 ) than does the northern South China Sea (100.09 gC 
m~ 2 a -1 ) 48 , directly supporting a higher biomass of zooplankton 
and fish 19,4 '. Finally, the lower species diversity in the East China 
Sea (4,167 species), compared with that in the South China Sea 
(5,613 species) 50,51 , likely reduces interspecific competition and 
allows for increased population sizes and carrying capacity. Given 
these factors, the present-day differences in effective population sizes 
of cutlassfish suggest a close correlation to the primary productivity 
differences in the East China Sea and northern South China Sea. 

The historical sizes of T. japonicus populations in the East China 
Sea and northern South China Sea during the period of original 
coalescence and first growth peaks are comparable (Fig. 4), suggest- 
ing similar carrying capacities for these two regions prior to 45 kya. 
This finding is supported by the geological history of the East China 
Sea and northern South China Sea during the Pleistocene, specifically 
the wide continental shelf subsidence 150 kya 18,52 (Supplementary 
information S2). Given their habitat dependence on the continental 
shelP 4 , the narrow continental shelves of the East China Sea and 
northern South China Sea prior to 150 kya would have provided 
limited habitat and supported similar population sizes of cutlassfish 
observed here. 

The similar population sizes of T. japonicus observed in the East 
China Sea and northern South China Sea 45-150 kya (Fig. 4) are 
surprising given the formation of the wider East China Sea contin- 
ental shelf — 150 kya 52 . Given that present population size is closely 
correlated to primary productivity, this finding suggests that carrying 
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Figure 5 | History of monsoons and population dynamics of T. nanhaiensis and T. japonicus in the northern South China Sea over the past 165 kya 
(2.03% per myr). Top panel: East Asian monsoon changes recovered from aeolian quartz flux (EQF) and fluvial quartz flux (FQF) of Lake Biwa 55 . Bottom 
panel: BSPs of T. nanhaiensis and T. japonicus in the northern South China Sea. Abbreviations (EG and GB) were defined for figure 4. 



capacities and levels of primary productivity between the East China 
Sea and northern South China Sea during this period were compar- 
able. In terms of primary productivity, temperate rivers are thought 
to play more important roles than tropical rivers by contributing 
terrigenous organic matter and nutrients into these marginal seas 
(Supplementary information S3). In this context, the Changjiang 
contributes over 90% of river-derived runoff and particles into the 
East China Sea 53 and has likely had a large impact on the historical 
primary productivity of the East China Sea 19 . However, the geologic 
history of the Changjiang development during the late Pleistocene is 
unclear (Supplementary information S4). 

If the Changjiang has drained into the East China Sea since the late 
Pleistocene, differences between the T. japonicus population sizes in 
the East China Sea and northern South China Sea during the last 
interglaciation (70-140 kya) or interstadial period (30-60 kya) 
should be similar to those of the present, given the increased primary 
productivity of the East China Sea 19 and the reduced primary pro- 
ductivity of the northern South China Sea 54 induced by strong sum- 
mer monsoons and weak winter monsoons 55,56 (Fig. 5). However, this 
model is not supported by our data. 

The slight glacial bottleneck observed in T. japonicus from the East 
China Sea is also indicative of the Changjiang history. During the 
LGM (22.5-30 kya), the sea level fell by 150 m in the East China 
Sea 57,58 , greatly reducing the sea surface area and exposing much of 
the continental shelf (Fig. 1). If the Changjiang had flowed into the 
East China Sea prior to the LGM, freshwater and terrigenous matter 
flux would have greatly declined due to the reduced annual runoff 



from the lower Changjiang drainage by intense dry and cold winter 
monsoons (Fig. 5). With less nutrient input from the Changjiang, 
primary productivity in the East China Sea would have greatly 
dropped during the LGM. Subsequently, a severe population bottle- 
neck for T. japonicus relative to earlier population sizes would have 
been expected due to losses of primary productivity and habitat size 
in the East China Sea. However, this is not true based on the observed 
slight bottleneck during the LGM (Fig. 4), suggesting a limited 
decrease in the primary productivity of the East China Sea during 
the LGM relative to that prior to the LGM and supporting the idea 
that the Changjiang did not historically drain into the East China Sea 
until the postglaciation and Holocene periods. 

The recent abrupt population growth (over a 10-fold increase) 
observed in the East China Sea population of T. japonicus during 
0-22.5 kya [0-15 kya for 2.7% per myr, 0-30 kya for 1.35% per myr] 
(Fig. 4) lends support to the hypothesis that the drainage of the 
Changjiang shifted to the East China Sea during the postglacia- 
tion-Holocene and increased primary productivity. Evidence from 
sediments of the Changjiang Delta also supports this hypothesis 59 61 . 
The geologic history of the Changjiang Delta can only be traced back 
to the late Pleistocene (MIS 3) 59 or Holocene 60,61 . Therefore, runoff 
from the Changjiang could not have contributed a significant 
amount of organic matter to the East China Sea until sea levels rose 
during the postglaciation or Holocene period 61 , suggesting that the 
earlier primary productivity of the East China Sea was similar to that 
of the northern South China Sea. Based on these patterns, it appears 
that the demographic history of T. japonicus in the East China Sea 
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Late Pleistocene climatic cycles 



Glacial low sea level and intense winter 
monsoon 




Interglacial high sea level and intense 
summer monsoon 










Low precipitation and terrigenous 
influx into the East Asian marginal seas 




High precipitation and increased land-source 
influx into the East Asian marginal seas 










northern SCS: Increased palaeo- 
productivity by strong upwelling 
current; 

ECS: Low palaeo-productivity 




northern SCS: Low palaeo-productivity due 
to weak upwelling current; 

ECS: Increased palaeo-productivity through 
the postglacial Changjiang discharge 



Population dynamic response of cutlassfish in the China seas: 

Northern SCS: LGM growth and postglacial decline; 
ECS: LGM decline and postglacial growth 



Figure 6 | Correlation of climatic changes with the palaeo-productivity and population dynamics of cutlassfish in the China seas. In the East China Sea, 
increased precipitation and discharge from the Changjiang during the intensified summer monsoons caused an increase in primary production. 
However, in the northern South China Sea, primary production was mainly influenced by nutrient advection caused by upwelling currents during winter 
monsoons. Thus, different population dynamic responses to climatic cycles were observed between the East China Sea and the northern South China Sea. 



closely reflects the historical changes in the drainage of the 
Changjiang into the East China Sea since the LGM. 

Compared with the glacial bottleneck and recent rapid population 
expansion revealed for the East China Sea population of T. japonicus, 
both the T. japonicus and T. nanhaiensis populations from the north- 
ern South China Sea exhibited a pattern of glacial growth and post- 
glacial-Holocene decline in population size (Figs. 4, 5). Contrasting 
demographic patterns have also been revealed for regional popula- 
tions of European Nassarius nitidus between two glacial refugia: an 
unaffected slow exponential growth in the Adriatic population dur- 
ing the LGM and a dramatic exponential growth after LGM in the 
Iberian Atlantic population 62 . The unique oceanographic and eco- 
logical conditions (e.g., strong interspecific competition due to high 
species diversity 50 , fishing pressures and changing primary produc- 
tivity) may all have had dramatic impacts on the population sizes in 
the northern South China Sea. 

However, the stress of increased competition appears not to have 
caused synchronous interspecific bottlenecks during the postglacia- 
tion-Holocene. While the northern South China Sea population of T. 
japonicus experienced continuous population expansion during 45- 
165 kya, the population of T. nanhaiensis steadily declined (Fig. 5, 
Supplementary information S5). Despite the increased catch produc- 
tion of the past few decades 27 , no similar population bottleneck has 
been detected for T. japonicus in the East China Sea (Fig. 4). 
Furthermore, cutlassfish are not an economically important species 
in the northern South China Sea, and their catch productions exhibit 
a fluctuating cycle of ca. 10 years 63 , suggesting that fishing pressure 
has also not played an important role in influencing population 
dynamics of cutlassfish. 

The close correlation between present-day effective population 
sizes and primary productivity suggests that the historic palaeo-pro- 
ductivity along the northern slope of the South China Sea may be 
responsible for the similar glacial growth and postglacial/Holocene 
genetic bottlenecks exhibited by T. japonicus and T. nanhaiensis. In 
the northern South China Sea, primary productivity was historically 
mainly influenced by fluctuating winter monsoons 54 (Supplementary 
information S6). The peak population size of T. japonicus and the 



rapid expansion of T. nanhaiensis 22.5-30 kya coincide with these 
palaeo-productivity increases of the LGM. The subsequent rising sea 
levels of the postglacial period and weakened winter monsoons could 
have caused a population decline in T. japonicus. The contraction of 
T. japonicus coupled with the rising sea level of the South China Sea 
would have temporarily provided new niche space for T. nanhaiensis, 
leading to the continued expansion of T. nanhaiensis 15-22.5 kya. 
However, the T. nanhaiensis population also declined, limited by 
decreasing palaeo-productivity during the postglaciation-Holocene 
(0-15 kya) 64 . In this context, the similar glacial population growth 
and postglacial genetic bottlenecks experienced by T. japonicus and 
T. nanhaiensis in the northern South China Sea appear to indicate 
general responses to the palaeo-productivity history of the South 
China Sea during the postglacial stage and Holocene. Similar con- 
cordant responses of interspecific genetic bottlenecks to range-size 
constraints caused by Mid-Quaternary glaciation (0.6-1.1 Ma) were 
also observed in two Patagonian fish (Galaxias maculatus and G. 
platei) 65 . 

In summary, the use of BSP analyses and a calibrated Cyt b 
molecular clock of 2.03% per myr revealed a close correlation 
between Pleistocene glacial cycles, primary production and popu- 
lation dynamics of cutlassfish in the China seas. The early synchron- 
ous population growth of T. japonicus in both the East China Sea and 
the northern South China Sea likely indicates similar oceanographic 
conditions prior to 45 kya. However, the recent intraspecific or inter- 
specific demographic differences of T. japonicus and T. nanhaiensis 
may reflect different evolutionary responses to differences in regional 
primary productivities induced by palaeo-hydrology and East Asian 
monsoon changes between the East China Sea and the northern 
South China Sea (Fig. 6, Supplementary information S7). 
Therefore, our study suggests the likelihood of genetic signs of 
Pleistocene climatic cycles and the postglacial development of the 
Changjiang Delta. A small number of similar studies regarding the 
dynamic response of different populations to regional environmental 
conditions have been conducted to date (e.g., European marine 
Gastropoda 62 , Patagonian freshwater fish 4 and humans 2 ). However, 
future comparative studies on a wider scale regarding the evolution- 
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ary dynamics of marine taxa from different ecoregions are needed 
because they may provide additional information into past palaeo- 
climatic and oceanographic conditions. 

Methods 

Sampling and data collection. Specimens of T. japonicus and T. nanhaiensis were 
sampled from 1 1 fishing ports along the coast of the northern South China Sea 
(sampling sites 1-11, Fig. 1 ) . Additional specimens of T. japonicus were sampled from 
14 offshore stations in the East China Sea and one coastal fishing port in Xiamen, 
Fujian (sampling sites 12-26, Fig. 1). Muscle tissue samples from tail clips ( — 10 cm) 
were preserved in 95% ethanol for genetic analyses. Total genomic DNA was 
extracted from each individual using a standard phenol/chloroform approach as 
outlined in Sambrook et al. 66 . Approximately 10-30 ng of DNA was utilised as a 
template to amplify the full-length mitochondrial Cyt b gene (ca. 1200 bp) via 
polymerase chain reaction (PCR). Amplifications were conducted in 25 uL volumes 
containing 5 uL of template DNA, 1 X PCR reaction buffer, 2 mM MgCl 2 , 300 nM 
each of primers CytbF and CytbR 67 , 200 uM dNTPs, 0.5 unit of Taq polymerase and 
ddH 2 0. Reactions were conducted under the following PCR conditions: an initial 
denaturation at 94°C for 3 min, followed by 35 cycles of 94°C for 30 s, 59 C C for 60 s, 
72 "C for 70 s, and a final extension of 72 C for 10 min. PCR products were separated 
on a 1.5% agarose gel and purified with the Gel Extraction Mini Kit (Watson 
BioTechnologies, Shanghai, China). Purified products were sequenced in both 
directions using the primers pair CytbF and CytbR on an ABI Prism 3730 automatic 
sequencer. 

Data analyses. Genetic diversity and haplotype networks. Cyt b gene sequences were 
aligned using Clustal_X 68 under default parameters, and the quality of the alignments 
was checked through the correct translation at the amino acid level. To confirm 
species identification, sequenced haplotypes were compared with T. japonicus 
(DQ364150) and T. nanhaiensis (JX477078) sequences from GenBank. 

To visualise relationships among Cyt b haplotypes, intraspecific median- joining 
networks 69 were constructed for each species using NETWORK version 4.6.1.0 
(Fluxus Engineering, 2012). To avoid hyperlinks and complex reticulations between 
haplotypes, the weight of hypervariable sites (number of mutations > 12) were set to 
0, and the weights of other variable sites were proportionally adjusted according to 
their mutation frequencies (e.g., the weights of minimum mutation events (number of 
mutations = 1) were increased to 90, and the weight of sites with 10 mutations was 
down weighted to 9). 

To reveal the regional responses of T. japonicus to different oceanographic con- 
ditions, the T. japonicus data set was divided into two geographic groups corres- 
ponding to two main marine ecosystems/ecoregions: the northern South China Sea 
and the East China Sea. Though these seas connect through the Taiwan Strait, they 
possess different biogeochemical properties and species assemblages 23,70 and can be 
divided near the southern tip of Taiwan Island and Dongshan of Fujian province 
(dashed line in Fig. 1). Nucleotide (tt) and haplotype (h) diversity estimates were 
calculated for each species, sampling site and geographic group according to the 
methods of Nei 71 using ARLEQUIN 72 (version 3.5). 

Molecular rate calibration. DNA calibration methods based on fossil data and ancient 
geologic events (i.e., > 1 Mya) have been shown to be unreliable and to underestimate 
evolutionary rates at the population level 40,41,73 . However, several recent geological 
and climatic events may serve as possible calibration points for the evolution of 
cutlassfish in the China seas: ( 1 ) the formation of the wide continental shelf of the East 
China Sea 150 kya, when coastal East China subsided during the last uplift of the 
Tibetan Plateau 18 ; (2) the flooding of the newly formed East China Sea shelf during sea 
level highstands of the Last Interglaciation (70-140 kya) 74 ; and (3) the postglaciation- 
Holocene (0-22 kya) flooding of the East China Sea shelf as sea levels rose again after 
the Last Glacial exposure. Given the dependence of cutlassfish on continental shelf for 
habitat 24 , the newly formed wide East China Sea shelf and its subsequent flooding 70- 
140 kya likely prompted the earliest period of population expansion for cutlassfish in 
the East China Sea. Thus, an approximate molecular clock for Trichiurus was inferred 
from the population expansion parameter t (x = 2ut; where t measures time in 
generations, t measures time in units of l/2u generations, and u is the mutation rate of 
the entire region of DNA under study) 29 ; where t = 70-140 kyr and t (t = 4.289) is 
estimated from the mismatch distribution analyses in ARLEQUIN. The mutation rate 
parameter, u, was estimated from the equation u = mu, where m is the fragment 
length (m — 1127 for Cyt b) and u is the within- lineage mutation rate per nucleotide. 

Demographic history. Historical demographic patterns for each species and for the 
two geographic groups of T. japonicus (northern South China Sea and East China Sea) 
were investigated in ARLEQUIN using the mismatch distribution of pairwise nuc- 
leotide differences between haplotypes. Populations that have experienced recent 
demographic expansions are expected to exhibit a unimodal mismatch distribution, 
whereas a multi-modal distribution is expected for populations in genetic equilib- 
rium 29,75 . Furthermore, two neutrality tests, Tajima's D 76 and Fu's F s 77 were used to 
infer past population dynamics. Both methods provide information regarding 
demographic history in the absence of selection, with significant negative or positive 
values generally suggesting population expansions or bottlenecks 76 78 , respectively. 

To compare complex Skyline-Plot models against constant demographic models, 
extended Bayesian Skyline Plots (EBSPs) 79 were run to estimate the number of 
changes present in each data set. The demographic histories of T. japonicus and T. 



nanhaiensis were further investigated using a coalescent -based BSP model with a 
Markov Chain Monte Carlo (MCMC) sampling procedure in BEAST v 1 .7.4 80 under a 
strict molecular clock, 2.03% (1.35%-2.7%) per myr. Three independent runs were 
combined using LogCombiner and visualised in TRACER 81 . The MCMC analysis was 
run for over 6 X 10 8 generations to ensure a high effective sample size (ESS: >200), of 
which the first 10% was discarded as burn-in. The evolutionary model of each data set 
was estimated using Modeltest 3.7 82 . 
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